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The main purpose of this paper is to provide an asymptotically 
optimal test. The proposed statistic is of Neyman-Pearson-type when 
the parameters are estimated with a particular kind of estimators. It 
is shown that the proposed estimators enable us to achieve this end. 
Two particular cases, AR(1) and ARCH models were studied and the 
asymptotic power function was derived. 

1. Introduction. Local asymptotic normality (LAN) for the log likeli- 
hood ratio was studied for a several classes of nonlinear time series model, 
from a LAN the contiguity property follows, for more details the interested 
reader may refer to [2], [12], and [4]. Applying the contiguity property, we 
construct a statistic for testing a null hypothesis Hq against the alternative 

(n) 

hypothesis H\ , often a various classical test statistics depends on the cen- 
tral sequence which appears in the expression of the log likelihood ratio, in 
the case when the parameter of the time series model is known we obtain 
good properties of the test, precisely, the optimality, see for instance [9, The- 
orem 3]. However, in a general case, particularly in practice, the parameter 
is unspecified, in the expression of the estimate central sequence appears an 
additional term which is non degenerate asymptotically. The latter, alters 
the power function of the constructed test. 

In order to solve this very problem, and on a basis of an estimator of the 
unknown parameter, we introduce and define another estimator which does 
not effects asymptotically the power function of the test, more precisely the 
additional term is absorbed. The principle of this construction is to modify 
one of the component of the first estimator in order to avoid the additional 
term, the details of this method are expanded further in the section 2. 
The main purpose of this paper is to investigate the problem of testing two 
hypothesis corresponding to a stochastic model which is described in the 
following way. Let {(Yi,Xi)} be a sequence of stationary and ergodic ran- 
dom vectors with finite second moment such that for all z £ Z, where Yi 
is a univariate random variable and A,- is a cZ-variate random vector. We 
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consider the class of stochastic models 

(1.1) Yi = T(Z l ) + V(Z i )e u 

where, for given non negative integers q and s, the random vectors Zi is 
equal to (li-i, ^-2, ■ ■ ■ , Yi-si Xi, . . . , Xj_ g ), the ej's are centered i.i.d. 
random variables with unit variance and density function /(■), such that for 
each i G Z, a is independent of the filtration T% = a(Zj,j < i), the real- 
valued functions T(-) and V(-) are assumed to be unknown. We consider 
the problem of testing whether the bivariate vector of functions (T(-), V(-)) 
belongs to a given class of parametric functions or not. More precisely, let 

m = { (m( P , -),a(e, •)) , if/, ey & e x x e 2 } , 

Gi x 9 2 C R e x RP, 6i / 0, 9 2 / 0, where for all set A, A denotes the 
interior of the set A and the script " ' " denotes the transpose, I and p are 
two positive integers, and each one of the two functions m(p, •) and a(9, ■) 
has a known form such that a (9, •) > 0. For a sample of size n, we derive a 
test of 

(1.2) H : [(T(-),V{-)) G .M] against ^ : [(T(-),V{-)) £ M] . 
It is easy to see that the null hypothesis Hq is equivalent to 

(1.3) H :[(T(-),V(-)} = (m(po,-),<T(0„,O), 
while the alternative hypothesis H\ is equivalent to 

Hf.KT^Vi-)} / (m(p ,-),^o,-)), 
for some (p' ,9' Q )' G 61 x 6 2 . 

In the sequel, our study will be focused on the following alternative hypothe- 
ses. For all integers n > 1, the alternative hypothesis .ff| n ^ is defined by the 
following equation 

(lA)Hi n) :[(T(.),V(-)] = (m(p ,-)+ n- 1 2G(-),a(9 ,-)+ n-sSi-)), 

where G(-) and S(-) are two specified real functions. The situation is different 
in the case when the used statistic is the Neyman- Pearson test which is based 
on the logdikelihood ratio A n defined as follows 

(1.5) A n = logf A) = y>g(<7n,i), 

\jn,oj ~, , 
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where / n ,o(0 an d fn(') denote the probability densities of the random vec- 
tor (Y\, . . . ,Y n ) corresponding to the null hypothesis and the alternative 
hypothesis, respectively. 

The use of the Neyman-Pearson statistics needs to resort to the following 
conditions: 

Under the hypothesis Hq, there exists a random variable V n such that 

V n AAA(0,T 2 ), 

where — > denotes the convergence in distribution and some constant r > 
depending on the parameter = (p' ,0' Y, such that 

(1.6) A„ = V n (<Ao)-^y^+op(l). 

The equality (1.6) is a modified version of the LAN given by [9, Theorem 1]. 
We mention that there exist other versions of the LAN, we may refer to [11], 
[8], and the references therein. On the basis of the LAN, an efficient test of 
linearity based on Neyman-Pearson-type statistics was obtained in a class 
of nonlinear time series models contiguous to a first-order autoregressive 
process AR(1) and its asymptotic power function is derived (see, [9, Theorem 
1 and Theorem 3]). The expression of the obtained test depends on the central 
sequence V n (</>o) which itself depends on the parameter ^. In a general case 
the parameter 4>q is unspecified, so, in order to estimate it, we introduce, 
under some assumptions, an estimate preserving, asymptotically, the power 
on Neyman-Pearson test when we replace, in the expression of the statistics, 
the parameter 4>q by an appropriate estimator, <j) n . Say, this estimator will be 
constructed on the tangent space with the direction of the partial derivatives 
of the central sequences in <j) n , where 4> n is a -^/n-consistent estimator of 4>q . 
In the sequel, 4> n will be called a modified estimate (M.E.). 
This paper is organized as follows: Section 2 describes the methodology used 
to construct the M.E. In Section 3, we give the asymptotic properties of the 
proposed estimate. In Section 4, we conduct a simulation in order to evaluate 
the power of the proposed test. All mathematical developments are relegated 
to the Section 5. 

2. Estimation with modifying one component . Consider the prob- 
lem of testing the two hypothesis Hq against which are given in (1.3) 
and (1.4) respectively and corresponding to the stochastic model (1.1). We 
assume that the LAN (1.6) of the model (1.1) is established, for example 
refer to [9]. 
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Let (j) n = (p' n , 9' n )' a i/n-consistent estimate of the parameter 4> = (p' , 9' )' , 
where 

Pn = \Pn,l, ■ ■ ■ , Pn,e), 8' n = \ ^n,l, ■ ■ ■ ,0n,p) , 



po = (pi,...,pe) and 9 ' = (6q 



Our purpose is to construct another estimate <f>' n of the parameter (p' Q ;9' o y, 
such that the following fundamental equality is fulfilled 

(2.1) Vni^n) - V n (0 n ) = D n , 

where D n is a specified bounded random function. In the sequel, the func- 
tions 

(p, •) — > m(p, •) and (9, •) — > a{9, ■) are assumed to be twice differentiable. 
Our goal, is to find an estimate cj> n satisfying (2.1) pertaining to the tangent 
space r„, such that, for [X' ,Y')' el'x MP, the following equation holds 

r n :V n ((x,Y))-v n (0 n ) = dv' n (k)-((x -p n )',(Y -§ n )y, 

where 

'dVn^n) dV n (4>n) 9V n (4>n) ^nQ>r 



dVnM 



dpi dpi 89i 89 p 



and the script " • " denotes the inner product. 

With the connection with the equality (2.1), the new estimate is then given 
by imposing that the value (X',Y')' satisfied the following identity 

(2.2) D n = dV n ^ n )'.((X-p n )\(Y -9 n )'^j . 

Clearly, the equation (2.2) has I + p unknown values, so it has an infinity 
of solutions, after modification of the J n -th component of the first estimate 
p n , we shall propose an element in tangent space T n which satisfies the 

equality (2.2). We obtain then a new estimate <f>' n = 0n = (p'ni^'n)' °f 
the unknown parameter (J)q, where 

Prl = \PnAi ■ ■ ■ iPnl 



and such that: for s € {1, ... ,£}, p n , s = p n , s if s ^ j n and p n j n / p n>jn . 

The use of the notation <j) n explains that we obtain the new estimate ^ 
of the parameter <j)Q when we change in the expression of the estimate (, 
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the j n component with respect to the first estimate p n corresponding to the 
step n of the estimation. It follows from the equality (2.1) combined with 
the constraint (2.2) that 

8=1 ° Ps t=l ° Ut 

( 2 - 3 ) = — a" [Pn,jn ~ Pn,jJ- 

U Pjn 

By imposing the following condition 
(2.4) dVniM + 

and with the use of the equality (2.2) combined with (2.4), we deduce that 

(2 - 5) Pn ^ = ^jt) +pnJn ' 

In summary, we define the modified estimate by 

4>'n = 4>n^ = (j>n,l, Pn,j n -l, Pn,j n , Pn,j n +l, ■ ■ ■ , Pn,£, &n,l, ■ ■ ■ , #n,p) ■ 

With a same reasoning as the previous case and after modifying the k n - 
th component with respect to the second estimate, we shall define a new 
estimate 

such that for t G {1, . . . ,p} 

On,t = @n,t ift^k n and 9 n ,k„ / n ,k n - 

we obtain 

(2.6) vM M )-Vn{k) = ^^-(e n , kn -L, kn ). 

Under the following condition 

(2.7) + 0, 
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it follows from the equality (2.2) combined with (2.7), that 

( 2 - 8 ) Qn,k n = Jx \ +9n,k n - 

OVn{<Pn) 

In summary, we obtain the modified estimate 

0'n = 4>n' kn ' 1 = \Pn,l, Pn,l, #n,l, • • • , &n,k n -l > @n,k„ , #n,fc„+l , ■ ■ ■ , #n,p) ■ 

The estimate (f>n'^ (respectively, 4>n' kn ^) is called a modified estimate in 
j n -th component with respect to the first estimate (respectively, in k n -th 
component with respect to second estimate), we denote this estimate by 
(M.E.). 

Remark 2.1. For each step n of the estimation corresponding a value 
of the position j n or k n of the component where the estimate was modified. 

3. Properties of the (M.E.). 

3.1. Consistency. Throughout, <p n is a -y/n-consistent estimate of the un- 
known parameter (f>Q. The conditions (2.4) and (2.7) are not sufficient to get 
the consistency of the modified estimate (M.E.). In order to get its consis- 
tency, we need to resort to one of the following additional conditions . 



(C.l) 

1 8V n 



\fn 9p jn 

(C.2) 

1 dV n (d 



p 

— > c\ as n — > oo, 



p 

— > C2 as n — > oo, 



fn d9 kn 

where c\ and C2 are two constantes, such that c\ ^ and C2 7^ 0. 

Our first result concerning the consistency of the proposed estimate is sum- 
marized in the following proposition. 

Proposition 3.1. Under (2.4) and (C.l) ((2.7) and (C.2), respectively), 

the estimate 4>n^ n ^ (4>n ,kn \ respectively) is a \fn- consistent estimator of the 
unknown parameter <f>Q. 

In practice, it is not easy to verify the condition (C.l) (respectively, (C.2)), 
in the case when the unknown parameter cfto is univariate, a sufficient con- 
dition will be stated in Lemma (3.1), in this case, we need the following 
assumption: 



ASYMPTOTICALLY OPTIMAL TESTS 



7 



(C.3) : For all real sequence ij]n) n >\ with values in the interval [0, 1], we 
have: 

—^Vn(Vn(pO + (1 - Vn)4>n)) = Op(l), 



n 

where V n is a second derivative of V n - 

Remark 3.1. In a problem of testing the two hypothesis Hq against 
Hi , and when the error e, 's are centered i.i.d. and eo — > A/"(0, 1), a large 
classe of time series model satisfied the condition (C.3), for instance, we 
cite the nonlinear time series contiguous to AR(1) processes , the details 
are expanded further later in the proofs of the Propositions (3.3) and (3.4) ■ 

Now, we may state the sufficient condition which implies assumptions (C.l) 
corresponding to the case when the parameter of the time series model is 
univariate. 

Lemma 3.1. Let <p n be a \Jn- consistent estimate of the parameter <f)Q. 
Let c\ be a constant, such that c\ ^ 0, then we have: 

(i) Under (C.3), if -^=V n (</> ) — > c\, as n ^ oo, then \/A > 0, 
1 



P 



n 



Vn{4>n) ~ Cl 



> A — > 0, as n — >■ oo. 



3.2. Absorbtion of the error. Consequently, with the modified estimate 
and in the case when the error between two central sequences is bounded, 
it is possible to absorb this error, this result is stated and proved in the 
following proposition. 

Proposition 3.2. Let 4>' n be an estimate (\fn consitency) of the param- 
eter (p',9'y. We assume that there exists a known bounded function D n , 
such that 

(3.1) V n (4> n ) = V n (<t>0)-D n + Op(l). 

Then, there exists an estimate (ft n of (p',0')' such that 

V n (<Pn) = V n (<M + Op(l). 

Remark 3.2. The equality (3.1) gives the link between the estimated 
central sequences V n {4> n ) and the central sequence V n {<f>o). Sometimes it is 
not easy to establish the form of the function D n , in the next section, we 
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propose, under some assumptions, how to specify this function in the two 
cases, i.e., the case when the problem of testing the linearity and nonlin- 
earity of the s-th order and the case of time series model with conditional 
heteroscedasticity respectively corresponding to the equalities 

Yi = PqY^x + aG(Y(i - 1)) + and 
Yi = po^-i + aG(Y(i - 1)) + y/1 + pB{Y{i - 1)) e u 



respectively, where Y(i — 1) = ^l^__i,li_2, . . . , li_ s J, a and ft are real pa- 
rameters and the ei's are centered i.i.d. random variables with unit variance 
and density function /(•). 

Throughout, we assume that the function /(•) is positive with a third 
derivative, we denote by /(■), /(•) and f^(-) the first, the second and the 
third derivative respectively. For all i£M, let 

M,(*) = m 

According to the notation (1.5), we suppose that the three following condi- 
tions are satisfied : 

• (L.l): maxi<j< n \g n>i - 1| = o P (l), 

• (L.2): there exists a positive constante r 2 such that 
E?=l(<fcM-l) 2 = r 2 + p(l), 

• (L.3): there exists a.-F n mesurable V n satisfying Y^i=i(9n,i — 1) = V n + 

o P (l), where V n Af(0, r 2 ). 

Conditions (L.l), (L.2) and (L.3) imply under Hq the local asymptotic nor- 
mality LAN corresponding to the equality (1.6), for more details see ([9, 
Theorem 1]). This last theorem is the fundamental tool used later to aim to 
establish the LAN for the considering models. 

3.3. Link between central sequences in nonlinear time series contiguous 
to AR(1) processes. Consider the s-th order (nonlinear) time series 

(3.2) Y i = p Y i _ 1 + aG(Y(i-l)) + e i , \ Po \<l. 

In this case and with the comparison to the equality (1.1), we have 

Z i = Yi , T(Zi) = Po Y i _ l + aG(Y(i-l)) and V(Z l ) = l. 

In the sequel, it will be assumed that the model is a stationary and ergodic 
time series with finite second moment. We consider the problem of testing 
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(n) 

the null hypothesis Hq : a = against the alternative hypothesis H\ : a 
n~2 5 with the comparison to (1.3) and (1.4), we have 

m(p ,^-i), <t(0 o , Fi-i)) = (po^-1,1), X = {m(p,-),p G 61} , 
Z{ = (Y i _ 1 ,...,Y i _ s ) and S(-) = 0. 



Note that this problem of testing is equivalent to test the linearity of the 
s-th AR(1) time series model when (a = 0) against the nonlinearity of the 
s-th AR(1) time series model when (a = n~5). 

Throughout, the scripts " || • \\e+ p " , " || • \\e " and " || • || p " denote the euclidian 
norms in and MP respectively. It will be assumed that the conditions 

(A.l) and (A. 2) are satisfied, where 

• (A.l): There exists positive constants 77 and c such that for all u with 
\\u\\e +p > i], G(u) < c\\u\\ i+p . 

• (A. 2): for a location family {/(ej — c), — 00 < c < —00}, there exist a 
square integrable functions $1, ^2 and a constant 5 such that for all 
€i and |c| < 5, such that : 



d k f(e l - c) 



f(ei) dc k 



<^>k(ei), for k = 1,2. 



Under the conditions (A.l) and (A. 2) the LAN of the time series model (3.2) 
was established in ([9, Theorem 2]), the proposed test T n is the Neyman- 
Pearson statistic which is given by the following equality 

Tn = H^y- Z(a) }' r2 = E ( M /( e o)) E ( G2 ( y (°)))' 

and Z(a) is the (1 — a)-quantile of a standard normal distribution $(•). In 
this case, the central sequence is given by the following equality 



Vn(po) 



1 n 

-_^M / (e i )G'(y(z-l)), where r 2 = B(M 2 f (e ))B(G 2 (Y(0))), 
v n i= i 



and such that under Hq, V n (po) — > A/"(0, r 2 ). The asymptotic power of the 
test is derived and equal to 1 — &(Z(a) — r 2 ), recall that when po is known, 
this test is asymptotically optimal, for more details see [9, Theorem 3]. 
Our aim is to specify the form of the function D n which is defined in (3.1), the 
parameter po is estimated by the -^/n-consistent estimator p n and the residual 
ei is estimated by €i iTl = Yi — Yi_\p n . We have the following statement: 
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Proposition 3.3. Assume that the conditions (Al) and (A.2) hold and 
Ci 's are centered i.i.d. and cq — > N (0, 1). We have 



(3.3) 
where 
(3.4) 
(3.5) 



V(p n ) = V n (po) - D n + P (1), 



Pn 



Vn 



+ p n and c\ 



-c X s/n{p n - po) 



-E 



Y G(Y(0)) 



Remark 3.3. 



The use of the ergodicity of the model imposes to 



require the condition E 



< oo, therefore we choose the func- 



tion G(-) in order to get this condition. For instance, we shall choose 
G(Y(i - 1)) = where a ^ 0. 



• With this choice of the function G, the condition (Al) remains sat- 
isfied, in fact, we can remark that \G{u)\ < 2\a\, then for all u with 
\\u\\e+ P > T) we have G(u) < 2a x \\u\\ t+p x •p^j- < ^ x ||tt||^ +p , 
therefore, we shall choose c = — . 

j > n 

3.4. An extension to ARCH processes. Consider the following time se- 
ries model with conditional heteroscedasticity 



(3.6)y 4 = po^-i + a G(Y(i - 1)) + y/l + PB(Y{i - 1)) e u i £ Z. 

It is assumed that the model (3.6) is ergodic and stationary. It will be 
assumed that the conditions (B.l), (B.2) and (B.3) are satisfied, where 

• (B.l): The fourth order moment of the stationary distributions of (3.6) 
exists. 

• (B.2): There exists a positive constants rj and c such that for all u with 
\\u\\t +p > rj, B(u) < c||u||! +p . 

• (B.3): for a location family {b~ 1 f( Sj ^), - oo < a < -oo, b > 0}, 
there exists a square integrable function (p(-), and a strictly positive 
real <j, where ? > max(|a|, \b — 1|), such that, 



d 2 b~ l f(^) 



f(ei) do) db k 



where j and k are two positive integers such that j + k = 2. 
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We consider the problem of testing the null hypothesis Hq against the al- 

(n) 

ternative hypothesis H^ such that 

H : m(p, Zi) = p Yi^i and cr(0 o , •) = 1, 
H[ n) : m{p, Zi) = po^i-i + n~^G(Y(i - 1)) and a(6 , Z t ) = \Jl + rT^B{Y{i 

Remark that Hq, fl^ correspond to a = f3 = (linearity of (3.6)) and 
a = (3 = n~2 (non linearity of (3.6)) with the comparison to the equality 
(1.1), we have 



Zi = Yi, T(Z i )=p Y i _ 1 + aG(Y{i-l)) and V(Z t ) = y/l + f3B(Y(i - 1)). 
Note that when n is large, we have 



i 



a(0 , Zi) = sjl + n-^B(Y(i - 1)) ~ 1 + ^B(Y(i - 1)) = 1 + n~35(y(* - 1)). 

Under the conditions (A.l), (B.l), (B.2), and (B.3), the LAN was established 
in [9, Theorem 4], an efficient test is obtained and its power function is 
derived. In this case, the central sequence is given by the following equality 

1 f n n 1 

V n ( P0 ) = --j= r£ M f^i)G(Y(i - 1)) + £(1 + eiM f (ei))B(Y(i - 1)) j , 

such that under Hq, 

Vn(po) AAA(0,T 2 ), 

where 

r 2 = I E (G(Y(0)) 2 + ^"^ E (B(Y(0)) 2 + I±E (G(Y(0))B(Y(0)) 

where Ij = EfejAfJ(eo)) and j = 0, 1, 2. 
The proposed test is then given by 

(3.7) T. - /{^^W}- 

By the subsisting po by its y'ra-consistent estimator p n in the expression of 
the central sequence, we shall state the following proposition: 
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Proposition 3.4. Suppose that the conditions (A.l), (B.l), (B.2) and 
(B.3) hold and 's are centered i.i.d. and eo — > J\ (0, 1). We have 

(3.8) V{p n ) = V n (po)-D n + op(l), 
where 

(3.9) D n = -cx<Jn(p n - p ), 

(3.10) p n = J^L_ + p n a nd c x = -m\y o G(Y(0)) 

VnVPn) 1 

3.5. Optimality of the proposed test. Throughout, T n and f are the statis- 
tics test and the constant respectively obtained with the subsisting of the 
unspecified parameter (j>o by its modified estimate (j> n in the expression of 
the test (3.7) and the constant r appearing in the expression of the log like- 
lihood ratio (1.6) respectively. 

(n) 

We assume in the problem of testing the two hypothesis Hq against 
that the LAN of the the model (1.1) is established, in order to prove the op- 
timality of the proposed test. To this end, we need the following assumption 



(E.l) There exists a -y/n-estimate (j) n of the unknown parameter <j)Q and a 
random bounded function D n , such that 

Vn(<t>n) = V„(0o)- D n + P {\). 

It is now obvious from the previous definitions that we can state the following 
theorem: 



Theorem 3.1. Under LAN and the conditions (2.4) (respectively, (2.7)), 
(C.l) ((C.2), respectively) and (E.l) the asymptotic power of T n under 
is equal to to 

1 - $(Z(a) - f 2 ). 
Furthermore, T n is asymptotically optimal. 

We shall now apply this last theorem in order to conduct simulations 
corresponding to the representation of the derived asymptotic power func- 
tion. The concerned model is the Nonlinear time series contiguous to AR(1) 
processes with an extension to ARCH processes. 
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4. Simulations. In this section, we assume that e^'s are centered i.i.d. 
and e —> Af(0, 1), in this case, we have E(ej) = 0, E(e 2 ) = 1, and 
E(ef ) = 3. We treat the case when the unknown parameter <pQ = po £ ©i C 
K, under i?o, the considering time series model can also rewritten 

(4.1) Yi = poYi-i + ei where \p \ < 1. 

4.1. Nonlinear time series contiguous to AR(1) processes. To evaluate 
the performance of our estimator, we provide simulations with comment in 
this section. In the case when the parameter po is known, the test T n is opti- 
mal and its power is asymptotically equal to 1— 3>(Z(a)— r 2 ), for more details 
see [9, Theorem 3]. In a general case, when the parameter pq is unspecified, 
firstly, we estimate it with the least square estimates p n = l v Z — > sec ~ 

ondly, with the use of the (M.E.) under the conditions (2.4) and (C.l), the 
modified estimate p n exists and remains -y/n-consistent, making use of (2.5) 
in connection with the Proposition (3.3) it follows: 

(4.2) p n = ■ + p n = ^r-—, + Pn, 

with the substitution of the parameter pQ by its estimator p n in (3.7), we 
obtain the following statistics test 



Vn(Pn) 

r(pn) 



>Z{a)\ where f 2 = E(M / 2 (e ,„))E(G 2 (y )), 

and e ,n = Y - Y_ip n . 



It follows from Theorem (3.1) that T n is optimal with an asymptotic power 
function equal to 1 — <3?(Z(a) — r 2 {p n )). 

We choose the function G like this G : [x\, X2, x s , x s +i, x s+ 2, x s - 



^ where a ^ 0. 

In our simulations, the true value of the parameter po is fixed at 0.1 and 
the sample sizes are fixed at n = 30, 40, 80 and 400, for a level a = 0.05, the 
power relative for each test estimated upon m = 1000 replicates, we repre- 
sent simultaneously the power test with a true parameter po, the empirical 
power test which is obtained with the replacing the true value pq by its 
estimate (M.E.) p n corresponding to the equality (4.2), and the empirical 
power test which is obtained with the subsisting the true value po by its 
least square estimator LSE p n (an estimator with no correction), we remark 
that, the two representations with the true value and the modified estimate 
M.E. are close for large n. 




4.2. ARCH processes. With the substitution of the parameter po by its 
modified estimate p n , in (3.7), we obtain the following test 

such that 

f 2 = 7 , TO E (G(Y(0)) 2 + {h ' n ~ 1} E (B(Y(0)) 2 + Z 1|B E (G(Y(0))B(Y(0)) , 

4„ = E(4 in M|(e ,„)), J =0,1, 2, and e ,n = 31) - Y-lPn- 

In our simulations, the true value of the parameter po is fixed at 0.1 and 
the sample sizes are fixed at n = 30,40,80 and 200, for a level a = 0.05, 
the power relative for each test estimated upon m = 1000 replicates. We 

choose the functions G and B like this G = B : yxi,X2, • ■ -,x s ,x s+ i,x s+ 2, ■ ■ 
;x s+g j — > ^ where a / 0. 

We represent simultaneously the power test with a true parameter po and 
the empirical power test which is obtained with the subsisting the true value 
po by its estimate (M.E.) p n corresponding to the equality (4.2),we represent 
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simultaneously the power test with a true parameter po, the empirical power 
test which is obtained with the subsisting the true value po by its estimate 
(M.E.) p n corresponding to the equality (4.2), and the empirical power test 
which is obtained with the subsisting the true value po by its least square 
estimator LSE p n (estimator with no correction), we remark that, when n 
is large, we have a similar conclusion as the previous case . 



AR(1)-ARCH,n=30 AR(1)-ARCH,n=40 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



AR(1)-ARCH,n=80 AR(1)-ARCH,n=200 




Remark 4.1. We mention that the limiting distributions appearing in 
Proposition (3.3) and Proposition (3.4) depend on the unknown quantity 
b n = (f>n — Po)> *.e., in practice po is not specified, in general. To circumvent 
this difficulty, we use the Efron's Bootstrap in order to evaluate b n , more 
precisely, the interested reader may refer to the following references : [6] for 
the description of the Bootstrap methods, [1], [10] for the Bootstrap methods 
in AR(1) time series models and [7] for the ARCH models. 
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5. Proof of the results. 

Proof of the Proposition 3. 1 . Consider the following fundamental decom- 
position: 

(5.i) (^' jn) y = (kr + (o Jn y , 

where 0' jn = (Oj n ,i)' j +p y such that O jnji = when i ^ j n , 

and °jn,j n = Pn,j n - Pn,j n - 

p 

Firstly, we have 4> n — > secondly we can deduce from (2.5) that: 
(5-2) O jnJn = ^ T = ^D n 



dV„(4>n) \Jn _J_ dV„((f>n) 

9pj n y/n dp jn 



1 P 

Since D n is bounded, we can remark that -^D n — > 0, from (CI), there 
exists some constante c\ 7^ 0, such that dVn(0u) Cl; f rom (2.4) and 
since the function x — > - is continuous on M — {0}, it follows that the 

random variable J, ,, , — > — , then the couple ( -4= D n ; J, ,- , ) 

n/s 8 Pjn a Pjn 

converges in probability to the couple ^0 ; ^j-J , since the function (x, y) — > 
xy is continuous on 1 x 1, it result from (5.2), that the random variable 
°jn,j n £ = 0, therefore 

(5.3) (V = (0, ... 0, O in)in , ... 0)' A (0, ... 0, 0, .. . 0)'. 

Consider again the equality (5.1), since the function (x,y) — > x + y is contin- 
uous on W i+P xW i+p , it results from (5.3) that 4>n converges in probability 
to cf>o as 

n — > 00. Notice that the last previous convergences in probability follow im- 
mediately with the use of the continuous mapping theorem, for more details, 
see [3] or [13]. By following the same previous reasoning, we shall prove the 
consistency of the estimate (ftn' kn \ Note that (j)n^ n ^ is -^/n-consistent esti- 
mate of the parameter fa and 

-<h) = Op(l), 

where Op(l) is bounded in probability in In fact, it follows from (5.1) 

that 

(5.4) V^ Un) - <h) = Mk - <h) + V^o jn = o P {\) + V^o jn . 
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Since -i/nO,- ,- = D n J, ,- , and using the condition (C.l), it results 

that 

y/nOj n = Op 1 (l), where Op 1 (l) is bounded in probability in M. 
We deduce that 

(5.5) V^(^ Un) - 0o) = P (1). 

Notice that with a similar argument and with changing (j)n , (C.l) and 
(2.4) by 4> ( n kn \ (C.2) and (2.7) respectively, we obtain 

(5.6) ~ <h) = Op(1). 

In order to prove Lemma 3.1, we need to stated the following classical lem- 
mas: 



Lemma 5.1. Let pQ)ie{i,...,n. be a sequence of a positive random vari- 
ables on the probability space (Q,P,P), (ai)ie{i,...,«} a sequence of a positive 
(strictly) reals such that Yl\=i 7^ = 1> then we have, for each e > 0, 

Lemma 5.2. Let (X n ) n >Q be a sequence of a random variables on the 
probability space (fl, J 7 , P), such that X n = Op(l), then X% = Op(l). 

Proof of the Lemma 5.1. Firstly, we remark that, Ve > 0, we have 



E*>.}cy{*>i}. 



In fact, we suppose there exists 



\j2 Xi > e } and u * U { Xi > q 1 } 



then for each iE{l,. ..,/}, we have Xi(uj) < which implies that 

J2Xi(cj) <e, 

i=i 

hence a contradiction. With the use of the u-additivity, we obtain 



P 



'±*>')s*fe{*>i})*±*{*>i) 
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Proof of the Lemma 3.1. In this case (f>o = po £ ®i C R, we denote by 
p n the -y/n-consistent estimator of po. 

Let ^4 > 0, from the triangle inequality combined with the Lemma (5.1), we 
obtain: 



P 



-\=Vn{f>n) ~ Cl 

n 



> A 



= P 
< P 



n 



Vn(Pn) 



n 



Vn(po) 



+ 



J? 



Vn(P0 



Cl 



> A 



n 



Vn(Pn) 



n 



VniPo) 



> 



A 



+ P 



n 



Vn(po) - Cl 



> 



.4 



Firstly we have 



(5.7) 



P 



Vn(Po) 



Cl 



> 



A 



as n — > oo, 



Secondly we have 
1 



(5.8) 
(5.9) 



V n {p n ) 



n 



Vn(P0 
1 



n 



Vn(Pn 



Vn(Pn) 



Pn ~ PO 



where p n is a point between po and p n , then there exists a sequence % with 
values in the interval [0, 1], such that p n = 7] n po + (1 — rj n )p n , this implies 
that 

\p~n — Po\ < (1 — Vn)\p n — Po\ < |Pn — Po\, this last inequality enable us to 
concluded that p n is -y/n-consistency estimator of po, it follows from (C.3) 
applied on the equality (5.9) that 



(5.10) P 



1 



V n {p n ) 



1 



71 



Vn(po) 



> 



A 



as n — > 0. 



Thus we obtain (i). 



Proof of Proposition 3.2. It suffices to choose under (2.4) and (Cl) the 



estimate 



or under (2.7) and (C2) the estimate 



t(2,*!n) 



In order to prove the Proposition (3.3), we need a following classical result. 

Lemma 5.3. Let (Q,J-,P) be a probability space, (X n )n>i is a sequence 
of real random variables on ft. If X n converges in probability to a constant 
c, then, there exists a sequence of random variable (Y n ) n , with X n = c + Y n , 
such that, Y n converges in probability to 0. 
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Proof of Lemma 5.3. For all A > 0, we have: 

P {\Y n \ > A) = P (\X n - c| > A) ->■ 0, as n -> oo. 

Proof of Lemma 5.2. For all e > 0, 3Mi > such that: 
sup Q ((P (|X Q | > Mi) ) < e, this implies that sup Q ({P {\X a \ 2 > Mf) ) < e, 
therefore with the choice of M = Mi, we obtain the result. 

Proof of Proposition 3.3. ej's are centered i.i.d. and eo — > AA(0, 1), mak- 
ing use of the results of [9, Theorem 2], we have 

V n (po) = -^=^M f (€i)G(Y(i - 1)). 

j=i 

The estimated central sequence is 

1 " 

Vn(Pn) = —7= E - !))■ 

vn i=l 

By Taylor expansion with order 2, we have : 

(5.11) V n (p n ) - Vn(po) = V n (p„)(p n - p ) + ^V„(p„)(p n - p ) 2 , 

where p~ n is a point between po an d p n and 



-1 n 

v n (p n ) = -=^y i _ 1 G(y(i-i)). 



Note that 

Rn = l;Vn(p~n)(Pn ~ Pof = ^=^=V n (p~n) (\/™(pn ~ Po)) • 

Since the estimator p n is -y/n-consistent and with the use of Lemma (5.2), it 
results that 

v^QSn-po)) =O p (1), 
from the assumption (C.3), it follows that 

Rn = op(l), 

finally we deduce that, 

(5.12) V n (p n ) - V„(p ) = Vn(Pn)(Pn ~ Po) + Op(l). 
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This implies that 

Vn(Pn) Vn(po) V n (p n ) 1 V n (p'n) r( , . fU 

-{Pn ~ PO) + Op(l) = —j= -^y/n(p n ~ P0) + Op{l), 



(5.13) 

where p' n is between p n and po, and V n is the second derivative of V n . From 
the assumption (C.3), we have 

1 Vn(Pn) m 



>n 

since the estimator p n is -^/n-consistent, it result that 

VniPn) Vn(po) m 



Jn x/n 
this implies that 

(5.14) ^M = ^)+ 0P (1). 

In Jn 



With the use of (5.14), the equality (5.12) can also rewritten 

Vn(Pn) - Vn(Po) = ~ 7^ \/^(Pn ~ Po) + Op(l), 



(5.15) = M^l^Pn - Po) + op(1). 

/n 



It follows from the assumption (C.l) combined with the ergodicity and the 
stationarity of the model that, the random variable -^V n (po) converges in 
probability to the constant c\, as n — > +oo, where 



d = — E 



Y G(Y(p)) 



therefore from the Lemma (5.3), there exists a random variable X n , X n — > 
such that 

V n (p ) = Ci + X n . 

n 



We deduce from the equality (5.15) and the -y/n-consistence of the estimator 
p n , that 

(5.16) Vn(Pn) ~ Vn(po) = C\Jn(p n - p ) + Op(l) = -D n + Op(l), 

where D n = —ciJn(p n — po). Recall that the second derivative V n is equal 
to 0, this implies that the assumption (C.3) is satisfied. 
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Proof of Proposition 3.4- The assumption (C.l) remains satisfied and 
the proof is similar as the proof of Proposition (3.3), in this case, for all 
p £ 0i, we have 

-1 n 

^ n i=l 

By a simple calculus and since the the function / is the density of the 
standard normal distribution, it is easy to prove that the quantity 2Mf(p) is 
bounded, therefore, there exists a positive constant w such that 2MAp) < w, 
then 

1 1 " 

\^=V n (p)\<w-Y,yli\B(Y(i-l))\. 

Jn n 

i=l 

With the choice B(Y(i - 1)) = ^ with a ^ 0, it results that 

' i— 1 

l l n 

|-7=V„(p)| <2«;|a|-^y i 2 „ 1 . 

By the use of the ergodicity of the model and since the model is with finite 
second moments, it follows that the random variable ^ X^=i ^i-i ~^ ^> 
where k is some constant, this implies that the condition (C.3) is straight- 
forward . 

Proof of the Theorem 3.1. From the conditions (2.4) ((2.7), respectively), 
(C.l) ((C.2), respectively), it results the existence and the -y/n-consistency of 
the modified estimate estimate 4> n corresponding to the equation (2.5) ((2.8), 
respectively). The combinaison of the condition {E\) and the Proposition 
(3.2) enable us to get under Hq the following equality 

V„(0„) = V„(0o) + o P (l). 

This last equation implies that with op(l), the estimate central and central 
sequences are equivalent, in the expression of the test (3.3), the replacing 
of the central sequence by the estimate central sequence has no effect. LAN 
implies the contiguity of the two hypothesis (see, [5, Corrolary 4.3]), by Le 
Cam third lemma's (see for instance, [8, Theorem 2]), under H[ n \ we have 

V n AAA(r 2 ,r 2 ). 

It follows from the convergence in probability of the estimate <j) n to </>o, the 
continuity of the function r : • — > r(-) and the application of the continuous 
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mapping theorem see, for instance ([13]) or [3], that asymptotically, the 
power of the test is not effected when we replace the unspecified parameter 
(j>o by it's estimate, 4> n , hence the optimality of the test. The power function 
of the test is asymptotically equal to 1 — <$>{Z{a) — T 2 (4> n )), the proof is 
similar as [9, Theorem 3]. 
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